Effective time-independent description of optical lattices with periodic driving 
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For a periodically driven quantum system an effective time-independent Hamiltonian is derived 
with an eigen-energy spectrum, which in the regime of large driving frequencies approximates the 
quasi-energies of the corresponding Floquet Hamiltonian. The effective Hamiltonian is evaluated 
for the case of optical lattice models in the tight-binding regime subjected to strong periodic driv- 
ing. Three scenarios are considered: a periodically shifted one-dimensional (ID) lattice, a two- 
dimensional (2D) square lattice with inversely phased temporal modulation of the well depths of 
adjacent lattice sites, and a 2D lattice subjected to an array of microscopic rotors commensurate 
with its plaquette structure. In case of the ID scenario the rescaling of the tunneling energy, pre- 
viously considered by Eckardt et al. in Phys. Rev. Lett. 95, 260404 (2005), is reproduced. The 
2D lattice with well depth modulation turns out as a generalization of the ID case. In the 2D case 
with staggered rotation, the expression previously found in the case of weak driving by Lim et al. 
in Phys. Rev. Lett. 100, 130402 (2008) is generalized, such that its interpretation in terms of an 
artificial staggered magnetic field can be extended into the regime of strong driving. 

PACS numbers: 03.75.-b, 03.75.Lm, 67.85.-d, 67.85. Hj 



I. INTRODUCTION 



Optical lattices are artificial crystalline structures of 
matter prepared by subjecting ultra-cold neutral atomic 
gases to spatially periodic light-shift potentials arising 
in the interference patterns of multiple laser beams [l[. 
Quantum degenerate atomic samples arranged in opti- 
cal lattices allow to study tailored quantum many-body 
lattice models in a well controlled experimental environ- 
ment [31 . While a wealth of lattice geometries are nat- 
urally available, a variety of entirely new configurations 
and tuning options arises if in addition periodic driving 
is applied. Periodic shaking of a lattice, for example, 
permits to tunc the effective tunneling strength (even 
to negative values) [3, Q, which was recently used to 
drive a quantum phase transition between a superfluid 
and a Mott insulator p . Suitably tailored periodic driv- 
ing schemes allow to implement new building blocks for 
simulating electronic matter, as for example the effect 
of the Lorcntz-force acting upon the electronic charge in 
magnetic fields p-Q- This extends the scope of optical 
lattice models to include intriguing aspects of electronic 
matter, as for example, quantum Hall physics. 

A common approach in the analysis of periodi- 
cally driven quantum systems is to search for a time- 
independent effective Hamiltonian with an energy spec- 
trum approximating the quasienergies of the Floquet 
Hamiltonian of the system |9l-llll| . The accomplishment 
of this task typically requires to restrict oneself to specific 
classes of driving operators. In this article, an effective 
Hamiltonian is derived for an arbitrary driving operator 
in the regime of large driving frequencies. This effective 
Hamiltonian is evaluated for various driven optical lat- 
tice models. The lattices are assuined to operate in the 
tight binding regime described by Hubbard Hamiltonians 
p^ [l3| with additional external modulation. Firstly, a 
periodically shifted one-dimensional (ID) lattice is con- 



sidered. In the regime of strong driving, in accordance 
with previous work by Eckardt et al. [3|, a renormal- 
ization of the hopping amplitude J is obtained, which 
permits to tune J even to negative values, a scenario re- 
alized in a recent experiment [J] . In non-bipartite lattice 
geometries the selective adjustment of negative hopping 
amplitudes along certain directions of tunneling allows 
to simulate effects of frustrated magnetism [ij . Tuning 
of the energy associated with tunneling is a generic op- 
tion in driven optical lattice models, not easily realizable 
in solid state lattices. Secondly, a two-dimensional (2D) 
square lattice with inversely phased temporal modulation 
of the well depth of adjacent lattice sites is considered, 
and a similar rescaling of the hopping amplitude J is 
found. As a third example, yet unexplored in the regime 
of strong driving, a 2D square lattice is considered, which 
is subjected to an array of microscopic rotors commensu- 
rate with its plaquette structure. In previous work it was 
shown that for weak driving this staggered rotation acts 
to implement the effect of a staggered magnetic field, ap- 
lying flux with alternating sign to adjacent plaquettes 
] . Here, it is shown that for strong driving the struc- 
ture of the effective Hamiltonian and thus its interpreta- 
tion in terms of a staggered magnetic field is preserved 
up to non-local tunneling terms, which describe negligi- 
ble hopping between distant lattice sites. Similarly as 
in the first example, a rescaling of the tunneling energy 
arises. Our general expression of the effective Hamilto- 
nian should prove useful for analyzing further cases of 
interest. 

The article is organized as follows: in Sec. [Hi a few 
relevant elements of Floquet theory are recalled and the 
connection between the Floquet Hamiltonian of a general 
periodically driven system and the corresponding time- 
independent effective Hamiltonian is established. The 
effective Hamiltonian is expanded into a series of nested 
commutators. The resulting general expression is applied 



to the periodically shifted ID optical lattice in Sec. lIII Al 
to the 2D optical lattice with temporal modulation of 
the well depths in Sec. IIIIB| and to the 2D lattice with 
staggered rotation in Sec. lIVI Finally, the article is closed 
with conclusions in Sec. [V] Some straight forward but 
technical calculations are deferred to the appendix. 



II. FLOQUET DESCRIPTION AND EFFECTIVE 
HAMILTONIAN 



According to Floquet's theorem an arbitrary Hamilto- 
nian H{t) with periodic time-dependence (-ff(i) = H{t + 
T)) operating in some Hilbert space H possesses a set 
of T-periodic (i.e., \un{t)) = \un{t + T))) Floquet states 
\un{t)) and a spectrum of quasi-energies En determined 
by the eigenvalue equ ation fo r the Floquet Hamiltonian 



3 ecmation loi 



S){t)\Un{t)) =En\Un{t)) 



(1) 



The states |w„(t))e~'^"*/'' form a complete set of so- 
lutions to the Schrodinger equation Sj{t)\^p) = 0. Re- 
call that for each \un{t)) and arbitrary integer m the 
state \un{t),m} = |w„(t)}e'™"* with fl = 2tt/T is itself 
a Floquet state with quasi-energy En^m = En + mMl. 
The solutions to Eq. ([T|) thus display a Brillouin zone- 
like structure with respect to the time-axis with En to 
be chosen within the first zone [—Ml/2,Ml/2\. The 
states \un{t),m) with _E„ G [—h£t/2^'h£t/2] form an or- 
thonormal basis in the composite Hilbert space "H ® Ht, 
where TLt is the space of T-periodic complex-valued 
functions. Hence, Sn^n'Sm.m' = ((w„(i),m|w„/(t), m'))T 
with ((0(t)|^(t)))T = j. J^ dt{(f>{t)\^{t)))T denoting the 
scalar product in T-L^T-Lt- Equation ([T|) may thus be con- 
sidered as an ei gen value problem in the composite Hilbert 

spaceT^^Hr Ui. 

For an arbitrary stationary orthonormal basis \n) of 
% and an arbitrary time-periodic Hermitian operator 
F = Ft with F{t) = F{t + T) one may define an or- 
thonormal basis of the composite Hilbert space T-L (8 T-Lt 
by \n{t),m) = UF,m{t)\n) with Up^mit) = 
and arbitrary integer m. Defining 



-iF(t)+imQ,t 



^ 



ul, 



xmitpF.n^-it) 



iJeff^ {f,f"')T, 



(2) 



(3) 



where {■ ■ ■)t denotes time averaging over the periode T, 
it is straight forward to verify 



{{n{t),m\Sj{t)\n'{t),m'))T^ 

Sm,m' {{n\Has\n') + mMl) 
+ (1 - 5m,n') (n|(e^("'-")^^*Jo(?'°>)T|n') 



(4) 



According to Eq. ^ the matrix elements of S^{t) in 
the basis \n{t),m) of the composite space T-L ® Ht 



display a block structure with diagonal (to ~ to') 
blocks {n\Hcfi\n') + mh^ energetically separated by 
multiples of Kl and off-diagonal (to ^ to') blocks 
(n|(e*^"' -m)otj^(__. i^^jj^/^ coupling different diagonal 
blocks. Following a well-known result of perturbation 
theory, if 



„i(m'-m)Oty|(0,0) 



)T||<^f^ 



(5) 



is satisfied for all to 7^ to' with \\A\\ = Max{|(n|^|77i)| : 
n, to}, the off-diagonal couplings can be neglected yield- 
ing the approximation 



((n(t),TO,|i3(f)|n'(f),TO'))T» 

5ni.m' {{n\Hcs\n') + mMl) . 



(6) 



Thus, within the range of validity of condition ^, and 
given that the energy levels resulting from blocks with 
different to do not mix, i.e., if 



HcftW < Ml 



(7) 



holds, equation ([6|) shows that the quasi-energy spectrum 
of H{t) within the first Brillouin zone (to = to' = 0) co- 
incides with the energy spectrum of the time-averaged 
effective Hamiltonian ifeff- Consequently, within the 
sub-space associated to the first energy band En.m=o 
the time-dependent Hamiltonian H{t) and the time- 
independent Hamiltonian H^s are equivalent. 

The crucial task in practical applications is to identify 
a suitable operator F compatible with the constraints 
imposed by conditions ([5]) and ([7]). A useful recipe in 
this respect is to decompose the Hamiltonian H{t) ~ 
i/< {t)+Hy (t) into a weakly driven part i/< (t), satisfying 
ll-f^<(^)ll ^ ^^i ^-iid a strongly driven part H^{t), and 
then to choose F in order to integrate out only Hy(t), 

i.e., hF{t) = L dsHy{s). A useful expansion of S)p' 
in terms of multiple commutators involving F and F' = 
dF/dt can be derived. Defining the multiple commutator 
between operators A and B of order n-|- 1 by the recursion 

[A,B]^j^^ = [A,[A,B]J and [A,B]g = B, the practical 
relations 

,-n+l 



dt " ^An + l)\ 

n— 



£t;7TI)iI'^''''1' 



(8) 



JF 



G. 






hold for arbitrary operators F and G [l^, [l^ • Inserting 
Eq. © with G = H{t) into Eq. (g]) for to = to' = yields 



(9) 



n=0 



„ n\ \ n + 1 



Henceforth, the decomposition H{t) = H^{t) + Hy{t) 
is applied and F{t) is chosen to satisfy hF'{t) — Hy{t) 



One may then rewrite Eq. ([9]) as 

^(0.0) ^ 



(10) 



v^ * 



^-^ nl 



-nF(i),i?<(t)]„ 



-[F(i),i?>W]„ 



In the following sections expression ([TU|) will be applied 
to several examples of driven optical lattices and condi- 
tions ([5]) and d?]) are inspected to determine the range of 
validity for approximating -ff (t) by -ffcff- 



III. DYNAMICAL CONTROL OF TUNNELING 
IN OPTICAL LATTICES 

A. Periodic shaking of ID optical lattice 

First, the periodical shaking of a ID optical lattice 
is considered, which allows to suppress tunneling and 
even simulate negative tunneling energies. This scenario 
has been previously investigated theoretically in Ref. Q 
and experimentally in Ref. [J] and thus permits a use- 
ful test of Eq. ([TU| . The Hamiltonian is written as 
ii{t^ = Ho + W{t). Within the tight binding regime the 
time-independent Hamiltonian is the ID Bose-Hubbard 
Hamiltonian Hq = —JT^ + Hint with the tunneling oper- 
ators T± = ^ / , cj,c^ ± cj^Ci, and the onsite interaction 

Hint = ■jZ^i/^'^C^i' ~ 1) [Il)[ll- The periodic driving 
operator reads W{t) = 2Q cos{V,t) with Q = xJ2u ^"-i^- 
Here, c,y denotes the bosonic anihilation operator at site 
;/, fii, is the corresponding particle number operator, and 
(i/, fi) indicates summation over pairs of nearest neighbor 
sites. The parameters J, U, and x quantify the tunnel- 
ing strength, the on-site repulsion energy per particle, 
and the modulation strength, respectively. The operator 
Q acts to introduce a constant gradient of the chemi- 
cal potential and thus a constant force x/^^i where d is 
the lattice constant. Hence, the driving term W{t) repre- 
sents a tilt of the lattice with harmonic time-dependence. 
Experimentally, W{t) is realized by periodically shifting 
the ID standing wave forming the optical potential and 
transforming to the co- moving frame of reference [2l| . 

It has been shown recently in Ref. [3| that for suffi- 
ciently high driving frequencies the driven system H{t) 
behaves similarly as the undriven system Hq, but with 
the tunneling matrix element J replaced by the effec- 
tive matrix element JJo(2x//iJ7), where Jq denotes the 
Bessel function of order zero. Notably, since Jq can take 
negative values, negative values of the effective tunneling 
strength should become possible, a prediction confirmed 
experimentally [J|. This result is readily reproduced by 
means of Eq. ^. Choosing F{t) = 2Q am{nt)/{hn) 
yields hF'{t) = Wit) and [F{t),W{t)] = 0. Eq. ^ 
with i7< (t) = Hq and i7> (t) = W{t) thus simplifies to 



f,T' 



ri=0 



1 /2isin 



m{nt) y 



Ml J 



[Q,Ho 



(11) 



Use of [Q,Hint] — and evaluation of the commutators 
\Q,T±] =xT^ and thus [Q,r±]„ = x"T(-i)"± leads to 



.(0,0) _ u 7 Y^ 1 f2ixsm{nt) 



Hint — J / ^ 



n=0 

oo 



(2n)! V fi^ 



2n 



T+ (12) 



^ (2n-f 1)! V hn 



2n+l 



^^ 1 ^ 2^xsin(m) \ 

'^^(2n-fl)'' — ' 



Applying the relations 



(e'2™-sin2"(z)) - ( l)"(27i)! 






, (13) 
(-l)™(2n-M)! 



2^"(ri -|- m)\(n — to)! 



(n-|- m -I- l)!(n — to)! 

(e'2™^sin2"+i(z)) ^ (e'(2™+i)^sin2"(z)) ==0, 

for integers n, ?Ti < and (n — ?7i)! = if n < ?7i and 
(...) denoting the average over I-k with respect to z, and 
making use of the power expansion of the 77i-th order 
Bessel function 



°° / -lyi 2n+n 



n=0 



22"+'^n!(TO + n)l 



(14) 



finally yields 



^i2(m'-m)ntfi(0,0)v _ 



(e-- ~"-'^T')t 



= -J J, 



S. 



2{m' —m) 

Hu 



2X I rj. 



m' ,rn -^int 5 



(15) 



/ i(2(m'-m)-|-l)Ot^(0,0)v j j I --->- 1 ti 

(e ^ ^ '^ ' S^F IT = J J2(m'-m) + l \ T7^ ] T- 



2X 



Ml 



and (setting m ^ m') 



2x 



HcS = —J Jo { TTT ) T+ + Hint ■ 



(16) 



Equation (|16p confirms the rescaling of the tunneling en- 
ergy derived in Ref. [3|. Finally, the range of applica- 
bility of HcS remains to be discussed via inspection of 
Eqs. ([5]) and ([T])- For increasing number of particles 
in the optical lattice, these conditions avoiding any cou- 
pling between different Floquet bands become increas- 
ingly hard to satisfy, because multiple excitations (e.g., 
multiple particle-hole excitations due to collisions) then 
may bridge increasing energy intervals. However, as has 
been pointed out in Ref.[22|, such higher order excita- 
tions yield only small couplings acting on time-scales too 
long to be relevant in typical experiments. Thus, it suf- 
fices to require that the single particle energy scaling fac- 
tors I J J™ {2x/fin) I and U in Eqs. ^ and Eqs. ^ do 
not exceed Ml. Because Jm{z) < 1 for arbitrary inte- 
gers TO, a sufficient condition for any value of the driving 
strength x ^^ J ^ ^^ and U <^ Ml. 



In the derivation of Eq. ([T6| a harmonic modulation 
W{t) oc cos(r2i) has been assumed, however, other time- 
dependences yield a similar renormalization of the hop- 
ping strength J. If, for example, cos(flt) is replaced 
by a rectangular modulation function f{t) = (—1)'' for 
— 7r/2 + vn < iot < 7r/2 + z^tt, where i> runs through all 
integers, the zero order Besscl-function Jo hi Eq. (|16p is 
to be replaced by the sine function. 



B. Well depth modulation in 2D optical lattice 

Next, a 2D square optical lattice is considered, com- 
posed of two sublattices as illustrated in Fig. [TJ The 
difference between the well depths of the A and the B 
sites is assumed to be modulated harmonically. Experi- 
mentally, this scenario arises if two standing light waves 
with wavelength A and parallely oriented linear polariza- 
tions are crossed by superimposing the two branches of 
a Michelson interferometer. Control of the optical path 
length difference of the interferometer allows to adjust 
the difference 9 between the time-phases of the two stand- 
ing waves. The resulting optical potential is V{x, y) = 
-Vq \e'^/^sui{kx) + e-^^/'^sm{ky)\^ [H with k = 27r/A. 
It is straight forward to implement the temporal modu- 
lation 6* = (1 + KCOs{Q,t)) 7r/2. For n <^1 this yields an 
overall potential V{x, y, t) = —Vq (sin^(fcx) -I- sin^(fcy)) — 
Vqtt sm{kx) sin(fcy) cos(rii) composed of a stationary 2D 
square lattice with mean well depth Vq and the desired 
harmonic modulation of the well depth difference of A 
and B sites. Within a tight binding description in terms 
of a driven Hubbard model a calculation similar to those 
in Ref. [1, [13, EM yields the Hamiltonian 



H{t) 



T± 



Hu 



-JT+ + Hi, 



VK(t) 



reAi^61,2,3,4 



u 



y^ fir (fir - 1) 



(17) 



W{t) = 2Qcos0t, Q^-f X! 



tj. '''i--|-ei 



Here, a^ and 6r+e„ denote the bosonic anihilation opera- 
tors at sites r and r -I- e,y of sublattices A and B, respec- 
tively. The corresponding particle number operators are 



7t. 



-^t 



rir = a^Gr and fir+ei = ^r+ei^r+ei- The four vectors e^, 
;/ = 1, 2, 3, 4, connect an A-site to its four nearest neigh- 
boring B-sites, according to Fig.[TJ In analogy to the pre- 
vious section, the parameters J and U quantify, respec- 
tively, the tunneling strength and the on-site repulsion 
energy per particle of the conventional Bose-Hubbard 
Hamiltonian for a 2D square lattice. The modulation 
strength is x = ^ ttkVq J dxdy \w{x, 2/)|^ cos(fcx) cos(fcy), 
where w{x, y) denotes the Wannier function of the lowest 
band. Although the definitions here refer to a 2D lattice. 




FIG. 1: (color online) Decomposition of 2D square lattice 
into two sublattices indicated by A (red) and B (blue). The 
vectors e^, v = 1,2,3,4, connect an A-site to its four nearest 
neighboring B-sites. The grey area denotes a A/2 x A/2-sized 
plaquette. 



a close formal analogy to the previously discussed ID case 
can be observed. In fact, with F = 2Q sm{^t) / {Hi^) 
and the definitions of Eq. (J17p the same commutation 
relations as those found in subsection IIII Al arc recov- 
ered, i.e., [F{t),W{t)] = 0, [Q,i/int] = and [Q,T±]^ = 
X"r'{-i)'«±- Consequently, equations ([T^ . p3|) and ([TC)) 
also hold for the 2D scenario of Eq. ((T7| . As in the ID 
case, one may thus use the modulation parameter x in 
order to suppress tunneling or even adjust negative val- 
ues of the tunneling strength. Arguments analogue to 
those used at the end of subsection IIII Al show that for 
any value of the driving strength x the time-independent 
effective description is justified if J, C/ <C hQ. 



IV. 2D OPTICAL LATTICE WITH 
STAGGERED ROTATION 



In this section the square optical lattice in Eq. ([T7|) of 
Sec. IIII Bl is considered, however with a geometry of the 
driving term, specifically designed to apply angular mo- 
mentum with alternating sign to neighboring plaquettes. 
It has been shown in Refs. [2J, [2^ that this can be ex- 
perimentally realized in a bichromatic optical lattice, pro- 
duced in an optical set-up comprising two nested Michel- 
son interferometers. This yields an optical potential 
V{x, y, t) = Vl{x, y)+V[{ix, y, t) consisting of the station- 
ary square lattice VL{x,y) = — Vo (sin^(fca;) -I- sin^(/i:y)) 
of Eq. (|17p and a temporal modulation VE.{x,y,t) = 
KVi,{x,y) coa{2S{x,y) — &t) with k adjustable within 
[0, 1] and 



S{x,y) 



_ _^( sin(fca;) — sm{ky) 1 
[ sm{kx) + sm{ky) J ' 



(18) 



which acts as an array of microscopic rotors, each cen- 
tered in an individual plaquette. The well depth Vq scales 
linearly with the overall intensity of the lattice beams 
while K is adjusted via the intensity ratio of the two fre- 
quency componets of the bichromatic lattice. A descrip- 
tion in terms of a driven Hubbard model leads to (cf. 
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FIG. 2: (color online) Plot of J and |^a/| versus the lattice 
well depth Vb for k = 0.4 and k — 1. All quantities are 
scaled to the recoil energy _Ercc = h^k^/2m of the photons 
providing the optical lattice potential. 4 J is approximated by 
the energy width of the lowest band. 



Ref. m 



H{t) = -JT++Hint + W{t), 



W{t) = ^AriVcos(0i)-6/Afsin(0i), 



^-E( 



Jlj- ^r+ei > 



(19) 



rel 



M 



^ {-ir+\a%+,,^+R.c.) 



reA, 1^=1-4 



The operators T± and _ffint arc the same as in 
Eq. p7|) . The modulation strength parameters are 
given by £,n = 2kVo J dxdy |w(x, j/)p cos(fca;) cos{ky) and 
^M = kVq J dxdy'w*{x+\/4:,y)[sm^{kx)—cos'^{ky)]w{x~ 
X/A,y). Since the 2D Wannier function of the lowest 
band of the square lattice Vi,{x,y) factorizes {w{x,y) = 
w{x)w{y)), one obtains the simplified expressions ^n = 
2kVo [J dx\w{x)\'^ cos{kx)] and S_m = KVoJdxw*{x + 
X / i) sin^ {kx)w {x — X/ 4). Because w{x) is peaked around 
X = 0^ ^N is on the order of 2kVo. In contrast, ^m 
exhibits the same order of magnitude as J because it in- 
volves a tunneling integral, which receives its main con- 
tributions from the small side lobes of the Wannier func- 
tion. This is confirmed by a band structure calculation 
including the first five bands as is shown in Fig. [2l where 
J and \^m\ are plotted versus the lattice well depth Vq 
for K = 0.4 and its maximum possible value k = 1. For 
K >« 0.4 the value of £,m can exceed that of J for in- 
creasing Vq, which would lead to negative overall hopping 
amplitudes of H{t) for certain fractions of the modula- 
tion cycle. This indicates that the description in terms 
of the lowest band Wannier function of the stationary 
lattice Vl becomes questionable then. 

For weak driving (^m <C fifl and £,n <C hil) a cor- 
responding time-averaged effective Hamiltonian H^e = 
-JT+ - i/\i [M,N] with K = ^N^M/ifi^) has been in- 




1 st order 
transition 



2nd order 
transition 



FIG. 3: (color online) Zero temperature phase diagram of the 
Hamiltonian -JT+ - iK\ \M,N] + Hint within the {J,K)- 
plane calculated according to Ref. [7| for fixed chemical po- 
tential IJ./U. Within a rectangular region (grey) around the 
origin of the {J ,K)-p\a,ne a Mott insulator is formed. It is sur- 
rounded by four superfluid phases, each with an order param- 
eter characterized by different values at the four lattice sites 
at the corners of a plaquette as indicated by the small asso- 
ciated pictograms. The width of the Mott region is A{fi/U), 
where A{z) = [z{2g - 1 - z) + g{g - 1)]/[4(1 + z)] for a filling 
fraction g £ {1,2, ...} and z £ [g — l,g]. Solid (dashed) black 
lines denote second (first) order phase boundaries. 



troduced in Ref. [7|. The interest in this scenario arises 
because this Hamiltonian mimics the action of a stag- 
gered magnetic field upon charged particles. The mag- 
netic field arises from the operator i [M,N], which pro- 
vides tunneling with imaginary hopping amplitudes 



i[M,N] 



2i 



E 



{(-l)''aj6r+e„-H.c.} . (20) 



Notably, ioi K Ri J the magnetic flux per plaquette can 
be on the order of a fundamental flux quantum, a regime 
accessible in solid state lattices only at the several 100 
Tesla level. It has been shown in Ref. [Tj that the Hamil- 
tonian — JT+ — iA'i [M, N] + Hint possesses a rich zero 
temperature phase diagram. Depending on the tunneling 
parameters J and K, four superfluid phases with distinct 
symmetries of the corresponding order parameters can 
be formed (see Fig. [3]). Unfortunately, since S,m < J and 
^7v/ft0 ^ 1, the staggered vortex superfluid predicted in 
Fig. [3] if I A' I > \J\ cannot be accessed in the weak driv- 
ing regime. Furthermore, since J is positive, the entire 
left half-plane of the phase diagram in Fig. |3] cannot be 
explored. This raises the question whether strong driv- 
ing permits to access an extended portion of the phase 
diagram or rather gives rise to entirely new physics. 

In the following, Eq. (|10p is used to explore the regime 
of strong driving for two different choices of the operator 
F{t). Firstly, the regime of strong driving with respect to 
^nN cos(0t) is explored, i.e. ^a? ~ Ml while J, [/, S^m ^ 



hft. The choice H<{t) = -~JT+ + Hint - 6f^sin(0i), 
F>(i) = £,NNcos{0t), and F = ^n N sm{0t) / hn yields 
hF' = H>(t), [F,i7>(<)] = and [F, iJi„t] = 0. With 
these settings Eq. pH)) becomes 



— '/t 

iO^^^i^int - J^-sin"(0t) 



n=0 



h0 



N 



[N,T+ 



(21) 



-sin"+i(0i) 

7? ' 



n=0 






W,M]^ 



K,„(2).Jn(z 



-=^ Z 




Using the relations [iV, T±] = -2T:p and [TV, [TV, M]] = 
4Af shown in the appendix, the multiple commutators in 
Eq. ^ are calculated to be [N,T±]^ = (-2)"T±(_i),. 
and [TV,M]„ = 2"-[i-(-i)"l/2[Ar,M][^_(_-^j„j/2. Insert- 
ing these into Eq. (PT|) and making use of Eqs. ([T^ and 
(full leads to 



(e' 



-JJ2 



i2mntp,(0,0) 



)t 



5m, ^int 



26v 



T+ - i-^MK2r, 



26v 



[N,M] , 



(22) 



and 



„i(2m+l)atf,(0,0)v _ 



^\ 



'h 



where Km{z),m G {0, 1, ...} is defined by 

B f 1 + f-lV" z™ \ 

/.„W.|(j„W-^V^5^)^ (23, 

In particular, for m=0 one obtains 

Since for z <^1 (corresponding to weak driving) Ji{z) « 
2^/2 = ^Ar/?i0, Eq. (pi)) reproduces the result H^s ~ 
-JT+ - K\ [M, N] found in Refs. 0,13. In the regime of 
strong driving the formal structure of the effective Hamil- 
tonian is maintained, however with rescaled tunneling en- 
ergies. The scaling functions Jo{z) and Ko{z) = —Ji{z) 
illustrated in Fig. 2] (solid lines) are bounded with max- 
imal values smaller than unity. They exhibit zero cross- 
ings, which permits to suppress Jeff = JJq{z) and A'off = 
KKq{z) or even adjust negative values. 

Upon use of the arguments discussed at the end of 
subsection IIIIAl equations (|22)) and (|24)) imply that 
the description in terms of iJcff is suitable if the values 
of the parameters J, C/, S,m are constrained by the rela- 
tions JJjn{z),U,£,MKm{z) <^ h@. Since Jm{z) < 1 and 



FIG. 4: (color online) Plot of scaling functions K2m{z) for 
771 = 0, 1, 2, 3, 4, 5, 6 and Jo{z). 



K2m+i{z) < 1/2 for all m = 0, 1,2, ..., it is sufficient to 
require J,U,£,M,£,MK2m{z) < h0. AUthough K2m{z) is 
not bounded, one may allow K2m{z) to be significantly 
larger than unity without affecting the necessary condi- 
tion ^MK2m{z) ^ fi.0, since S,m <SC ^0- This in turn per- 
mits values of z significantly exceeding unity, i.e., within 
the regime of strong driving. This is illustrated in Fig. |4j 
where K2m[z) is plotted for m = 0,1,2,3,4,5,6. One 
recognizes that, in the area around z sa 2, where |_ftro(2)| 
attains its maximum, the values of K2m{z) with m > Q 
remain well below one, i.e., the requirement to prevent 
transitions between different Floquet bands is satisfied. 
In the region 2 < z < 2.5 the effective hopping amplitude 
Jeff can be tuned close to zero, while Kq[z) maintains siz- 
able values. This is the regime of interest in experiments 
exploring the staggered vortex superfluid phase of Fig. |31 
where tunneling is dominated by the staggered magnetic 

flux (Jfcff > Jeff). 

Although, the rotor potential defined in the context of 
Eq. P^ constrains S,m to be on the order of or smaller 
than J with the consequence that S,m ^ ^0, it is nev- 
ertheless instructive to examine an alternative approach 
to an effective Hamiltonian, which does not a priori pre- 
suppose small values of ^m ■ This amounts to the choice 
iJ<(t) = -JT+ and fl">(t) = W{t). For simplicities sake, 
the interaction TJint is neglected here. It turns out prac- 
tical to rewrite the driving term as 



W{t) = Qte<0t 4_ Qg-»0t ^ 



(25) 



With F = i(Q'fe«s^* - Qe-'^^^)/hQ. and thus hF' = W 
one finds [F, W] = 2i [Q, Qtj /^jj^^ which inserted into 
Eq. (dni) yields 



^T' - -JT.h^F^T,i 



(26) 



2 x^ ^- 



— Y 



nn ^-^ n!(n-f2) 

n=0 ^ ' 



[^,[Q,Ql], 



Evaluation of iJcft requires to calculate the multiple com- 
mutators [F,r+]„ and [F, [Q,Q^]]„- 



Calculation of [F,T^]^: a calculation outlined 
the appendix yields the commutators [A/, T±\ 

(1 T 1) EreA,.= l-4(-l)''+' («K+2e„ - 6U>r+e 

While this implies that [M, T+] = 0, one notes 
that [A/, r_] exclusively comprises terms denoting non- 
local tunneling between distant lattice sites. Such 
terms should not significantly contribute and are hence- 
forth neglected, and thus [A/, T_] « is assumed. 
In addition, the commutator relation [N,T±\^ = 
(-2)"r±(_i)„, already applied below Eq. ^, holds. 
Thus, the summand proportional to Af in _F ~ 
[N^N sin(fit) — M^M cos{ilt)] /Ml can not contribute to 



[F,r. 



+J 



-K 



N 



nn 



sin" {nt)T, 



(-1)- 



(27) 



Calculation of [F,[Q,Q'^]] : since [Q, Q^] = 
— ^^N^AiiNjM], multiple commutators of the form 
[F, [N, A/]]„ have to be considered. A direct calcu- 
lation (cf. appendix) yields [N, [N, M]] = AM and 
[A4^, [A/,iV]] = 167V 4- Onl, where Onl denotes the 
non-local tunneling terms specified in Eq. (|A2p . which 
are henceforth neglected. Introducing the abbreviations 
S = ^[N, M] and L = \M yields the relations 



[N,L] = S 
[N,S] = AL 
[L,S] = -AL 



(28) 



Because F = p^L +pnN is linear in L and A^ with p^ = 
— ^^cos(ili), and pn = |^sin(r2i), the multiple com- 
mutator [F, S']„ is a sum of commutators [Fi, . . . [F„, S]] 
where F^,i> G {!,..., n} denote either of the opera- 
tors plL or pnN. According to the operator alge- 
bra of Eq. (pS)) . for even n all possible combinations 
[F„, . . . [Fi, 5]] are scalar multiples of S, while for odd 
n scalar multiples of A' or L can arise. Furthermore, 
only those contributions [F„, . . . [Fi, S]] are non-zero, for 
which, except for Fn in case of odd n, the F^^s occur in 
pairs, i.e, -P^i^-i = ^2i^- Summing over all contributions 
yields after some algebra 






{pI 



-pI) 



s 



(29) 



{pI 



P%) 



{pnL -~plN) 



where n 6 {0, 1, ...}. 

With the expressions of Eq. (p7|) and Eq. (|29|) we are 
now in the position to calculate H^ff by means of Eq. (|26p . 
Upon noting that the time-averages of the odd com- 
mutators yield zero contributions, i.e., ([i^, T+Jj^ i ]^)t = 
{[F, [Q, Q^]] 2„+i)t == 0, and that for arbitrary x, y e M 



((a;2sin2(fii) -|-y2cos2(m))")T 
'2{n- v)\ f2v 



(30) 



E 



-V) 
V 



2;\ 2{n-iy) /yx 1v 

2/ V2. 



one obtains 




^cff = 


-^cffT'+ - i ATcff 


^eff = 


«(«) 


A'off = 





(31) 



Here, Jq{z) denotes the Bessel function of zero order and 

(_4)" _, 



Ko{x,y) = J2J2 



^ '^ (n + l)(2n)! 



(32) 



2{n-v)\ f2v\ /a;\2("-i') fy^^2y 
n-v )\ V fvl) V2 



Note the symmetry relation KQ{x,y) = Kq{ii,x). 
Furthermore, Kq{x^x) ~ sin(2a;)/a; — sin {x)/x'^ and 
A"o(0, y) = Ji{2y)/y. As a consequence of the latter 
equation, if Kq is constrained to the y-axis, one recov- 
ers the case of small ^m ^ Ml described by Eq. (|24p . 
The scaling functions Jq{z) and KQ{x,y) illustrated in 
Fig. [5] are bounded with maximal values of unity arising 
at the origin. In contrast to Jos, K^ff is not bounded 
with respect to the driving parameters ^m and ^^v since 
K itself scales with £,m^n- In order to determine the 
range of ^m and ^n for which conditions ^ and ([7]) 
hold, (e*(™'-'")"*.f)^'°V can be readily calculated by 
means of Eq. ^^ and variants of Eq. ^U^. With the 
same relaxed variants of conditions ([5]) and ^ used in 
the previous discussion one finds that, either of the pa- 
rameters S,m/MI and £^n /Ml may exceed unity, as long 
as their product remains well below unity. To illustrate 
the area of permissible values for S,m and ^jv, the red line 
in Fig. [5jb) shows the upper boundary of the region de- 
fined by ^M^N < 0.1(fi.r2)^. When passing too far into the 
region above this line, one may encounter transitions be- 
tween different Floquet bands such that the description 
in terms of the effective Hamiltonian Hcff fails. Recall, 
that at two points in the calculations leading to Eq. ([HT]) 
non-local terms have been neglected, which describe tun- 
neling between distant lattice sites. These terms arise in 
quadratic or higher order with respect to S,M/fi0, such 
that in the calculation of Eq. ([M|) they were a priori ex- 
cluded by the assumption ^m/^0 ^ 1- Finally, if the 
Hubbard interaction Hint had been accounted for in the 
derivation of Eq. (|3T|) . the choice HF' = W would have 
led to a non-zero commutator [F, H-mt] with the conse- 
quence of non-local interaction terms also scaling with 
£,M/h0 to second or higher order. 



V. CONCLUSIONS 

In this article, a general expression for the effective 
Hamiltonian of a periodically driven quantum system has 



y 2 




FIG. 5: (color online) (a) Plot of scaling functions Jo{z) (red, 
solid), Ko{z,0) (blue, dashed), and Ko{z,z) (black, dotted), 
(b) Plot of scaling function Ko{x,y). The thick red line is 
given by xy — 0.2. 



been derived and applied to various optical lattice mod- 
els with external driving, which are readily accessible in 
experiments. Three scenarios have been considered: a 
periodically shifted ID lattice, a 2D square lattice with 
anticyclic temporal modulation of the well depths of ad- 
jacent lattice sites, and a 2D lattice, subjected to an ar- 
ray of microscopic rotations commensurate with its pla- 
quette structure. For the ID scenario, the rescaling of 
the tunnelirig strength previously found by Eckardt et 
al. in Ref. [3| is reproduced, showing that strong driv- 
ing permits suppression of tunneling and the adjustment 
of negative tunneling energies. Well depth modulation 
inversely phased for adjacent lattice sites lets us extend 
these effects to a 2D scenario. In the case of the 2D square 
lattice with staggered rotation, the expressionpreviously 
found for weak driving by Lim et al. in Refs. [3,|3l is gen- 
eralized. The effective Hamiltonian is found to maintain 
its formal structure independent of the driving strength, 
i.e., the interpretation that it simulates the effect of a 
staggered magnetic field upon charged particles, previ- 
ously found for weak driving, extends to the regime of 
strong driving. For all lattice scenarios one finds that 
the regime of strong driving differs from that of weak 
driving essentially by an energy rescaling. In case of the 
2D lattice with staggered rotation, this rescaling should 
permit to extend the experimentally accessible portion of 
the phase diagram. The general expression for the effec- 



tive Hamiltonian might also be useful in numerous cases 
of interest not considered here. 
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Appendix A 



The following calculations are carried out for bosons. 



ar,a!, 



br,bl, 



Sr,r', while all other possible 

commutators are zero. Similar calculations yield analog 
results for fermions. The following definitions are used 



R 



E 



s = 



E 

reAi/el,2,3,4 



arbl+e^ 



(Al) 



reAi/el,2,3,4 



E i-^r'-rk 



r^^r+e^ 



Na =J2 «i»r ^^B^J2 ^rbr 



re/ 



relE 



R±R^ ,M = S + S^,N = Na 



Nf 



Calculation of [M,T±\: 
of Eq. jJT]) and [R, S] 

lows [A/,r±] = [S\R] T [S'^i?]^ One obtains 
five kinds of contributions to [S't,ii;]: 1) hop- 
ping terms connecting nearest neighbor A-sites 



with the definitions 
[i?t,S't] = it fol- 

t 



(-l)''+lat6r+e,.,< 



Bfi r+e„ 



= i-iy+^alc 



with odd 



/i ; 2) hopping terms connecting near- 



est neighbor B-sites (— l)''+-'^ajfer+e„, Or^r+e 

— (— l)''+^6r^g br+e^ with odd ly + ^i; 3) hopping 
terms connecting next nearest neighbor A-sites 

{-lY+^albr+^,^,ar+2eA+o,] = (-l)"+'4ar+2e,. ; 

4) hopping terms connecting next nearest neighbor B- 



i-ir+^albr+e^,arbl_,^^ 



sites 

5) on-site terms 



i^iy+'bl^^^br 



yir+^a%+,^,a,bl+^ 

(—1)"^+^ ( fltflr — 5r_|_e_^ 6r+e„ ) • Only the terms of type 3) 

or 4) can yield non-zero contributions, if the sum r g A 
and ly £ 1,2,3,4 is carried out, leading to [S't,/?] = 

EreA,i/=i-4(-l)''^^ (4ar+2e,. ^ 6j6r+2e J ■ Observing 

that [5"+,^]^ = [S^,R], one finally obtains [Af,T±] == 

(1 T 1) EreA,.=l-4i-^r^' (4«r+2e. " b%+2..) . 



Calculation of [N, T± 
Eq. (|AT|) one calculates [Na,T± 



With the definitions of 

T[Na,R]^. 
— —R one 



[Na,R 



With [Na,R] = EreA,i.=i-4 4ar,ar&I+, 

obtains [Na,T±] ~ —7^. A similar calculation gives 

[Nb,T±] = T± and thus [N,T±] = -~2T^. 

Calculation of [N,[M,N]]: Using the definitions of 
Eq. lAil) one calculates [M,N] = 2{S-S^). With 

[NaAM.N]] = 2{[Na.S] + [Na,S]^^ and [Na,S] = 

: —S one obtains 



Z^reA,i/=l-4' 1) 



i^+l 



[A'^, [A/, A^]] = — 2A/. An analog calculation yields 
[Nb, [M, N]] = 2M and thus [A^, [M, N]] = -AM. 

Calculation of [M,[M,N]]: With [M,N] = 
2(S'-S't) one gets [Af, [A/,A^]] = 4[5t,S']. 
Five kinds of contributions to [5*^,5] arise: 1) 
hopping terms connecting nearest neighbor A- 

sites [(-l)''+latfor+e„,(-l)^+'ar+e„+e,&l:+e„ 

— ajor+e^+e with odd ly + fi; 2) hop- 
ping terms connecting nearest neighbor B- 

sites 



{-ir+^albr+,^,{-ir+^arbl_^. 



Or+e„"r+e„ 



with odd 



M ; 3) hopping 



terms connecting next nearest neighbor A-sites 

(-l)''+14fer+e„,(-l)"+'ar+2e„fol+e„] = ala,+2.^ ; 

4) hopping terms connecting next nearest neighbor B- 



sites 



^r-e„^r+e„ ; 



5) on-site terms {-lY+^a%+^^,{-lY+^a,,hl 

ajor — fer^g^fer+e„ • Suming up over all r G A and 
v, ^ e 1,2,3,4 in compliance with the constraint of odd 



/i for type 1) and 2) terms yields 



[M, [Af , N] 
IQN + 4 



(A2) 



E 



flr^r+e^+e^, — ^r^r+e„+e 






The first term 16A^ results from the summation of the on- 
site contributions in 5), while all other contributions of 
type l)-4) are collected in the sum of non-local tunneling 
terms. 
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